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ABSTRACT 

Stellar mass loss rates are an important input ingredient for stellar evolution 
models since they determine stellar evolution parameters such as stellar spin-down 
and increase in stellar luminosity through the lifetime of a star. Due to the lack of 
direct observations of stellar winds from Sun-like stars stellar X-ray luminosity and 
stellar level of activity are commonly used as a proxy for estimating stellar mass loss 
rates. However, such an intuitive activity — mass loss rate relation is not well defined. 
In this paper, I study the mass loss rate of the Sun as a function of its activity level. 

1 compare in-situ solar wind measurements with the solar activity level represented 
by the solar X-ray flux. I find no clear dependency of the solar mass flux on solar 
X-ray flux. Instead, the solar mass loss rate is scattered around an average value of 

2 • 10"^"' Mq yr~^. This independency of the mass loss rate on level of activity can 
be explained by the fact that the activity level is governed by the large modulations 
in the solar close magnetic flux, while the mass loss rate is governed by the rather 
constant open magnetic flux. I derive a simple expression for stellar mass loss rates 
as a function of the stellar ambient weak magnetic field, the stellar radius, the stellar 
escape velocity, and the average height of the Alfven surface. This expression predicts 
stellar mass loss rates of I0~^^ — 10""'^^ Mq yr^^ for Sun-like stars. 

Key words: stars: magnetic field - stars: solar-type - stars: winds, outflows. 



> 1 INTRODUCTION 

X 



|Parker| ( |1958[ ) introduced the concept of a radially uniform 
supersonic flow, continuously flowing from an iso-thermal 
hot corona as a result of a hydrodynamic thermal accel- 
eration. One can imagine space as a huge vacuum cleaner 
that sucks up solar material, while solar gravity opposes 
this force. The gas slowly accelerates sub-sonically until it 
reaches a critical point, where its thermal energy equals the 
gravitational energy of the Sun. Above this critical point, 
the flow becomes super-sonic and continues to accelerate to 
some asymptotic terminal speed, Moo. This idea was contro- 
versial at flrst, but have been widely accepted since the first 
sporadic ( Gringauz et al.|1960 1 and consistent ( Neugebauer 



& Snyder 



1966 \ in-situ measurements of the solar wind were 



made, proving the existence of such flow. Parker's model has 
served as the baseline to all models for the solar wind and 
stellar winds; the main difference between the models is the 
source of the energy accelerating the wind. These models in- 
clude non iso-thermal winds, poljrtropic winds, sound waves 
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driven winds, dust driven winds, and line-drive winds (see 
[Lamers fc Cassinelli|199^ with references therein). 

Parker's model is hydrodynamic: the flow is assumed to 
be parallel to the magnetic field, and interaction between 
the flow and the fields is neglected. However, stellar mag- 
netic fields play a significant role in the physics of stellar 
winds since they govern the transport of the ionized gas in 
the corona ( Lamers & Cassinclli 1999|) . Parker ( 1958 1 did 
not include the solar magnetic field in his wind model, but 
he discussed extensively the consequences of the coupling be- 
tween the wind and the field to conclude the structure of the 
magnetic field in the interplanetary space. |Weber fc Davis| 
(1967) have shown that a star can lose angular momentum 
due to the fact that open magnetic field lines in stellar coro- 
nae (along where the wind is flowing) are attached to the star 
at one end, but are free at the other end beyond the Alfven 
point, VA, at which the flow speed equals the Alfven speed, 
UA ~ B/^/^np. Here B is the magnetic field strength and 
p is the plasma mass density. Coronal plasma is co-rotating 
with the star, but it is also frozen in to the magnetic field. 
Therefore, a torque is applied on the star by an arm with 
a length of va- This "magnetic braking" explains how stars 
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spin down, where the angular momentum loss rate depends 
on the stellar mass loss rate. Therefore, stellar mass loss 
rates and angular momentum loss rates are crucial inputs 
for stellar evolution models, which then provide input to 
other models. 

Unfortunately, neither stellar winds nor stellar mass loss 
rates can be directly observed for low-mass stars (seelGiidell 



2007 and references therein). Some attempts to estimate 



stellar mass loss rates were done using chemical separation 
and Hg profiles ([Michaud fc Charland|1986|[Lanz fc Catala 



1992 



1980 



1996 



Bertin et al. |1995l, radio observations ([Abbott et al 



ICohen et al ||1982| [Holhs et al.|pM5| |Lim fc White 
|Gaidos et al.|2000l, and observations of x-ray emission 



due to charge exchange (Wargelin & Drake 20011. These 



calculations have estimated stellar mass loss rates between 
10"^'' - 10"^° Mq yr. 

Perhaps the most extensive estimation of stellar mass 
loss rates was done using measurements of Lya absorptions 
at the edge of stellar heliospheres, where the stellar wind 



collides with the Interstellar Medium (ISM) (Wood et al, 
2002[ |2005| |Wood"2006). The formation of a shock at this 
location causes the build up of a Hydrogen "wall" that can 
be measured in the Lya line. By estimating the ISM prop- 
erties, the stellar mass flux carried by the wind can be de- 
termined. Based on many observations of different stellar 



types. Wood et al. ( 2005 I derived a scaling law for stellar 



mass loss rate as a function of x-ray luminosity, which can 
also be considered as a function of stellar age and activity 
level. They found that young, active stars (such as the young 
Sun) have mass loss rate of about 10~^^ Mq Yr~^ , which 
is ~100 times higher than the approximated solar mass loss 
rate of 2 ■ 10"^'* Mq Yr~^, derived from typical solar wind 
parameters near Earth. 

The results in [Wood et al.| ( |2005[ ) are presented cau- 
tiously, and the assumptions made are clearly stated. How- 
ever, their intuitive scaling law, which states that the more 
active the star, the higher the mass flux, is commonly used in 
an ad-hoc manner, despite of the fact that it might be over- 
estimated. Other calculations for the young Sun revealed a 
mass loss rate which is only 10 times higher than the cur- 



rent value (Holzwarth & Jardine 2007 Cohen et al. 2009 



Sterenborg et al. 20111. In fact, the scaling law in ( [Wood 
et al.||2005[) breeiks down for the most active stars. In con- 



trast, Fisher et al. (19981 and Pevtsov et al. (20031 used 



solar data to derive a scaling law between X-ray spectral 
radiance and magnetic flux. They found a nearly linear re- 
lation that holds for many orders of magnitude. Therefore, 
it is unclear how stellar winds and thus mass flux correlate 
with stellar magnetic activity and observed X-ray flux. 

Unlike the sparse information about stellar winds, the 
solar wind has been widely observed both remotely and in- 
situ since the beginning of the space era. These observa- 
tions, together with extended monitoring of the solar mag- 
netic field, have revealed a much more complex physics, with 
the solar magnetic field governing the power, the acceler- 
ation, and the topology of the solar wind. Magnetic en- 
ergization of the wind is necessary to explain the discrep- 
ancy between the fastest wind predicted by Parker's model 
(about 600 km s~^) and the observed fast solar wind (about 
800 — 900 km s~^). The exact mechanism in which the mag- 
netic energy is converted to kinetic energy is not fully under- 
stood yet. Models to explain this energy conversion consider 



the effect of wave dumping due to electron/ion resonance fre- 
quencies, dumping of Alfven waves, magnetic energy turbu- 
lent cascade, and energy dissipation due to magnetic recon- 

nection (see Aschwandcn 2005; Cranmer||2010 for detailed 
^ I II 1 1 II . I 

description and references therein). These models are tightly 
related to the solar corona heating problem. Great progress 
is currently being made in solving the coronal heating and 
solar wind acceleration problems due to new observations 
from the HINODEQand SDCQ missions. 

Based on the long-term in-situ measurements of the so- 
lar wind taken at 0.3 AU (HELIOS), the outer solar system 
(Voyager), high heliographic latitudes (Ulysses), as well as 
near Earth (WIND, ACE), some non-trivial properties of 
the solar wind have been revealed. These properties are well 



summarized in McComas (20071 (with references therein). A 
most notable feature is the bi-modal structure of the solar 
wind, which seems to be composed of two different popu- 
lations. The first population is the so-called "fast wind". 
It originates from coronal holes - the large-scale regions of 
open field lines (polar regions during solar minimum) - and 
it is rather steady, with speeds of the order of 750 km s~^ or 
more. The second population is the so-called "slow wind". 
This wind is much more sporadic than the fast wind and 
it is assumed to originate from the vicinity of the large, 
low latitude helmet streamers during solar minimum. Dur- 
ing solar maximum, the slow wind is observed to originate 
from higher latitudes, as well as from the vicinity of ac- 
tive regions. The two populations are not distinguished by 
flow speed only, but also by densities, temperatures, compo- 
sition, abundances, distribution functions, and First Ioniza- 
tion Potential (FIP) levels. While the acceleration of the fast 
solar wind could, in principle, follow Parker's model, with 
the contribution of magnetic energy by one of the proposed 
mechanisms above, the slow solar wind introduces a greater 
challenge in explaining its acceleration mechanism. An ob- 
served inverse correlation between the flow speed and elec- 
tron temperature, and a strong ion temperature anisotropy 
with very hot minor species at the low corona are inconsis- 
tent with the theory of thermal acceleration. In general, it 
seems that both populations are at first confined in closed 
loops, then are later released along open field lines. The fast 
wind is initially confined near coronal holes, where loops 
are small, cool, less dense, and short-lived, while the slow 
wind is confined within the large, hot, dense loops, with the 
plasma spending more time in the loop before it is released 
( [Feldmari eraT] [19991 ). |Fisk et al.[ (\^99\ showed how the 
mass flux and Poynting flux at the flux-tube base together 
with the lifetime of the loops are sufficient to explain the 
fast solar wind and determining its final speed. 

The difference in final speed suggests that the two pop- 
ulations gain different amounts of energy through the ac- 
celeration process. This can be generally explained by the 
different energy per unit mass the populations have (due to 
their different densities). In addition, due to its fast release 
the fast wind has more time to gain electromagnetic energy 
compared to the slow wind. The frequent ejection of the fast 
wind also explains why it is quasi-steady compared to the 
sporadic slow wind. Alternatively, the two populations can 
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gain the same amount of energy at the coronal base but have 
different amount of energy losses in the acceleration process. 
Schwadron & McComas (2003) proposed that both popula- 



tions gain similar electromagnetic energy to have high speed 
values (~ 800 km s^^). The difference in the final kinetic 
energy results from the wind radiating energy in the form 
of heat conduction and the advection of thermal energy as 
it is accelerated, where these radiative terms mostly depend 
on the maximum temperature of the loops. The loops can 
be thought as a reservoir of energy that is being transferred 
to the wind once it is released along open field lines. Since 
the slow wind is released from large, long-lived, hot loops, 
it radiates much more then the fast wind, so that its fi- 
nal kinetic energy is smaller. They also showed that their 
scaling law and its dependency on the maximum loop tem- 
perature explains the observed anti-correlation between the 
electron temperature and final wind speed, the FIP effect 
(Schwadron et al. 19991, and the anti-correlation between 



the magnetic flux-tube expansion factor and the wind speed 
( Wangfc Sheeley|1990|) 



'Schwadron et al. 



^(12006), and Schwadron & McComas 

(|2008J further showed that the energy function presented in 



Schwadron & McComas ( 2003 1 can be modified to provide 



a linear correlation between the solar X-ray power and the 
open magnetic fiux at the base of the fiux-tube. They showed 
that solar wind data obey this linear relation, that it can be 
extended to many orders of magnitude, and that it is similar 
to the relation presented in |Pevtsov et al.] ( |2003[ ). Recently, 
Wang (20101 have estimated the mass and energy fiux at 



et al. 



the coronal base using solar wind data taken by ACE and 
Ulysses in combination with solar magnetic field (magne- 
togram) data. Based on conservation laws along a magnetic 
flux tube, he derived an identical formula as in 
The only difference is that while 



et al. 



( 2006 

( 2006 1 assumed a constant value for the wind speed 



Schwadron 



Schwadron 



far from the Sun, Wang (20101 has used actual data for 
this parameter. He found that most of the energy fiux at 
the coronal base is used to lift the wind from the gravita- 
tional well, i.e. provides the wind with the escape velocity. 
He also found that while the magnetic field, the mass fiux, 
and energy fiux vary a lot at the coronal base, these param- 
eters are rather constant at 1 AU where the variations are 
balanced by the geometrical expansion of the magnetic fiux 



tube. Similarly, Schwadron et al. ( 2006 1 explained why the 



solar X-ray luminosity vary by an order of magnitude, while 
the observed wind power is rather constant. This is due to 
the fact that the strong solar X-ray emissions come from 
the closed, hot loops and the wind power is dominated by 
the open, rather constant solar magnetic fiux (Wang et al. 
2000[ ) . Both [Schwadron et al.| ( |2006[ ) and [Wang, ( ,2010[ ) have 



suggested that the solar wind fiuxes are determined by the 
large-scale, dipolar component of the solar magnetic field, as 
the heliosphere indeed seems to be composed by two hemi- 
spheres with opposite polarity separated by a single current 



sheet, which is reversed through the solar cycle (Smith et al 
2001| [Jones et al.|2003| . 

In this paper, I argue that the X-ray fiux observed on 
stars (or stellar activity in general) cannot be a reliable 
proxy for its wind mass fiux and power. In Section [2] I in- 
vestigate the dependency of the solar wind mass loss rate 
on the solar X-ray fiux using solar wind data taken at dif- 
ferent periods of time and at different locations in the solar 



system. In Section [3l I discuss these results and derive a 
relation between stellar mass loss rate, and the stellar mag- 
netic and gravitational parameters. I conclude my findings 
in Section m 



2 OBSERVATIONS OF SOLAR WIND MASS 
LOSS RATE AND SOLAR X-RAY FLUX 

In this section, I study the relation between the solar wind 
mass loss rate and the solar X-ray fiux, which is the proxy 
commonly used to relate stellar mass loss rates and the level 
of stellar activity. I use solar wind in-situ measurements 
taken by WIND/ ACE (near 1 AU), and by Ulysses (at high 
heliographic latitudes between 3-5 AU), as well as measure- 
ments of solar soft X-ray (1-8 A) fiux taken by GOES 7 
and ^ during solar cycle 23 (1996-2006). To study longer- 
term trends, I use solar wind observations for the years 
1965-201Crl Voyager II observations of the solar wind at 
the outer solar system, and long-term sunspot dataj (ISIDC- 
|team||2006| . ACE, WIND, Voyager, and Ulysses data have 
all been obtained from the CDAWEB websit<£] 



2.1 Data Selection 

Here, we are interested in comparing the global, ambient 
solar mass loss rate with a particular solar activity level. 
There are three challenges in this context when using in-situ 
data. First, transient events, such as Coronal Mass Ejections 
(CMEs) and Corotating Interaction Regions (CIRs) can ap- 
pear in the data, while they do not represent the ambient 
solar wind. Second, solar fiares and short-term significant 
increases in the X-ray data also do not represent the ambi- 
ent X-ray fiux. Third, an in-situ measurement of the solar 
wind does not necessarily represent the global value of some 
parameter. In order to avoid these issues, I select the solar 
wind data used in this study as follows. 

I first concentrate on the data available for solar cycle 
23 (taken between 1996-2006). I use hourly averaged data 
to obtain a large amount of data points. I exclude however, 
any impulsive, short-lived spikes so that the statistical be- 
havior is_dominated_b^;_tliearnbient values. Using the ICME 
list by Cane & Richardson (2003;rj I exclude periods which 
are associated with Interplanetary Coronal Mass Ejections 
(ICMEs) measured at 1 AU from the data. For Ulysses data, 
there is no such list. Therefore, I remove data periods when 
both the number density is higher than 10 cm~^ and the 
magnetic field magnitude is higher than 10 nT. Even though 
this method is robust, it is most unlikely to have such con- 
ditions in the ambient solar wind observed by Ulysses. 

The GOES data is provided in units of Watts/m^ as 
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observed at 1 AU. The X-ray flux units commonly used in 
the astronomy hterature are ergs/ {cm? s), where the X-ray 
luminosity is divided by the stellar surface. In order to be 
consistent with these units, I converted the GOES flux in 
to cgs units, and also used the factor {lAU)'^ / R^ to obtain 
the X-ray flux units in the astronomical context. When using 
the X-ray flux as the global representation for the solar level 
of activity, solar flares and short-lived data peaks must be 
removed. In order to exclude solar flares from the GOES 
data, I remove data points with X-ray flux higher than 
45 ergs cm~^ s~^ (10~® W/m^ - the typical flux associated 
with flares is above this value, see Golub fc Pasachoff|2009 1. 
Typical non-flaring X-ray fluxes in the GOES data are of 
the order of 0.1 -30 ergs cm"^ s"^ (5 -10"^ -10"^ W/m'^), 
which are three orders of magnitude lower th an the total 
solar X-ray flux - F^ kj 2 ■ lO" ergs cm"'^ s"^ ( |Wood||2006 
Golub fc Pasachoff|[2009| . 



2.2 Results 

Figure [l] shows scatter plots of the solar mass loss rate as 
a function of the solar X-ray flux for all valid data points 
measured near Earth and by Ulysses between 1996-2006. It 
can be clearly seen that the best linear flts have flat slopes, 
and while the X-ray flux varies by an order of magnitude, the 
mass loss rate is uniformly concentrated around the value of 
M0 « 2 • 10"" Mq yr-^ for aU X-ray flux values. The 
scatter around this constant value is slightly larger in the 
Ulysses data, probably due to the fact that Ulysses have 
measured more heterogeneous wind (fast and slow), while 
the 1 AU data contains mostly slow wind measurements. 

The best way to determine the global solar mass loss 
rate from in-situ data, is to combine simultaneous measure- 
ments taken at different positions. In particular, it is useful 
to use the available multi-point solar wind data taken at 
the same time by spacecraft near Earth at low heliographic 
latitude, and by Ulysses at high heliographic latitude. For- 
tunately, Ulysses was at high latitudes during 1996 (solar 
minimum) and 2000 (solar maximum). This enables us to 
investigate the properties of the slow and fast solar wind 
separately, while giving a good estimate of the total solar 
mass loss rate at a particular time. 

Figure |2] shows scatter plot of the wind speed and num- 
ber density as a function of the X-ray flux, only for data 
points taken during the year 1996 and the year 2000. Data 
points for the two years are shown in different colors. Fig- 
ure [2] clearly shows that during 1996, Ulysses has measured 
pure fast solar wind, while during 2000 (solar maximum), it 
has measured slow solar wind, and the interplanetary space 
was fllled with slow wind (even at high latitudes). During 
solar maximum, the large number of strong active regions 
appear on the solar surface take over the large-scale weak 
background dipole field component. As a result, much of 
the fast wind is eliminated and large amount of the solar 
wind originates from non-coronal hole regions ( [Schrijver fc| 
De Rosa|20"03} [Cohen et al.|2009[ ). 

In the four upper panels of Figure IS] I show the solar 
mass loss rate and the solar wind kinetic energy as a func- 
tion of the solar X-ray flux for data taken near 1 AU and 
by Ulysses during 1996 and 2000 (similar display as in Fig- 
ure [2| . The X-ray flux in Figure Isl as well as Figure p] has 
a clear separation between solar minimum period, with val- 



ues below 10 ergs s~^ cm~^ , and solar maximum period, 
with higher values ranging between 10 — 45 ergs s~^ cm~^ . 
The solar mass loss rate and the solar wind kinetic en- 
ergy are concentrated around 2 — 4 • 10"^'* A/q yr~^ and 
5 • 10^^ ergs s~^, respectively. The solar wind observed dur- 
ing 1996 was very steady and concentrated around these 
values, while the solar wind observed during 2000 is more 
sporadic, with larger variations around the mean value. In- 
terestingly, the fast wind seems to be slightly more energetic 
than the slow wind, while there seems to be no signiflcant 
difference in their mass fluxes. To provide an additional ref- 
erence observation, I also plot the solar mass loss rate and 
the solar wind kinetic energy as a function of time based 
on data taken by Voyager II during 1996 and 2000. At that 
time. Voyager II was located at the outer part of the solar 
system. Here I simply look at the values of these parameters 
throughout a year of a quiet Sun and an active Sun, without 
a direct comparison with the X-ray flux. The Voyager data 
shows variations of the solar mass loss rate, which are simi- 
lar to those plotted in the upper panels of FigurelS] and with 
a similar spread around the same average value. The solar 
wind kinetic energy is rather steady, but is slightly lower 
than the energy observed at inner parts of the solar system, 
in particular at the beginning of the year 2000. This can be 
due to solar wind energy exchange with pickup-ions, CIRs, 
or other temporal variations of the solar wind at that dis- 
tance, which cannot be identifled by this time-series alone. 

A longer-term analysis of the solar wind properties as 
a function of X-ray flux can be obtain by relating the X-ray 
flux with the most trivial solar activity proxy - the number 
of sunspots. I use sunspot data for the period when X-ray 
data is available to obtain a simple scaling of the X-ray flux 
with sunspot number (seen in the upper-right panel of Fig- 
ure ffl. The upper-left panels of Figure B] show scatter plots 
of the monthly averaged solar mass loss rate as a function 
of both the sunspot number and the scaled X-ray flux (sim- 
ilar to Figure [l| Figure [2] and Figure |3| for data taken near 
Earth between 1965-2010. It can be seen that here the slope 
is even more flat, and that when averaging the solar mass 
loss rate over a month, there is even less scatter around the 
same value of 2 ■ lO"^** Mq yr~^ . The two bottom panels of 
Figure |4] show similar plots as in the bottom of Figure |3J for 
daily and monthly averaged Voyager II data taken between 
1980-2009. These plots also indicate an average solar mass 
loss rate of 2 • 10"" Mq yr'^ . 

|Wang| ( |2010| ) has shown that there is a nearly perfect 
correlation between the mass and energy fluxes at the base 
of the corona. He stated that this is an indication that the 
majority of the wind energy is used to lift the plasma from 
the Sun's gravitational well. Figure [5] shows scatter plot of 
the mass loss rate as a function of the solar wind kinetic 
energy rate, both on a logarithmic scale, for data taken near 
1 AU and by Ulysses during 1996 and 2000, as well as sim- 
ilar plot for long-term solar wind data taken between 1965- 
2010 (monthly averages), and Voyager II data taken between 
1980-2009. The slope for all plots ranges between 0.5-1, sug- 
gesting that indeed, a more massive parcel of solar wind 
plasma has more energy. The spread and dynamic range of 
the plots is clearly related to the steadiness of the wind. The 
fast wind, observed by Ulysses during solar minimum covers 
a very narrow range, while the range is wider for slow wind 
observed near 1 AU during solar minimum. The range gets 
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even wider for observation of the slow solar wind taken both 
by Ulysses and near 1 AU during solar maximum. Naturally, 
the long-term data taken at 1 AU shows similar trend. The 
widest range appears in the Voyager II data, where in fact it 
is composed of two populations. One is the solar wind and 
the other is the gas measured at the Heliosheath. 



B\rA) ^ 
4TTp{rA) ' 



(1) 



or alternatively 



B {rAj = ATVp(rA)Ue = — ^~' 



(2) 



3 DISCUSSION 

Based on the results presented in Section [2. 2[ the mass flux 
carried by the solar wind does not have an obvious depen- 
dency on stellar activity. As mentioned in Section [l] this can 
be explained by assuming that the solar mass loss rate to the 
solar wind is determined by the amount of open magnetic 
flux opened up to the heliosphere. In contrast, solar activ- 
ity level is determined by the magnetic flux maintained in 
the closed loops. The order of magnitude variations of the 
closed flux through the solar cycle are much greater than 
the factor of 2 variations of solar open flux ([Wang et al 



|2000[ ). The small variations in the open flux are due to the 
fact that during solar minimum, most of the open flux orig- 
inates from the large polar coronal holes, while during solar 
maximum, much of it originates from the vicinity of active 



regions ( Schrijver fc De Rosa|[2003[ ). The amount of open 
flux is also modified by the level of magnetic reconnection 
between closed and open flux ( |Fisk|2005l ). This effect is en- 
hanced during solar maximum due to the complicated field 
structure, where the reconnection is driven by differential 
rotation, heliospheric current sheet tilt angle from the solar 
rotation axis, and CMEs which carry closed flux and force 



it to reconnect with the heliospheric open flux (Fisk 1996 



Crooker et aL][2002l [Owens et ar][2008l [20TT| [Yeates et aL] 



20101. As with the solar open flux, the mass loss rate of the 
solar wind is not constant, but is scattered around an aver- 
age value of about 2 — 4 • lO"'^* Mq yr~^ . Variations in the 
mass flux do not depend on variations in solar activity, but 
they are mostly due to density variations at the wind source 
region. 

Suppose the mass loss rate of a sun-like star is deter- 
mined by the amount of open stellar magnetic flux. Then 
we can assume that the wind mass flux is determined at the 
Alfven surface, where the stellar wind speed is equal to the 
Alfven speed, and all magnetic field lines are open. In addi- 
tion, the mass flux carried by the stellar wind needs to escape 
the gravitational well of the star, i.e. reach the stellar escape 
velocity, Ue = •v/2GM*/i?i, where G is the gravitational 
constant, and M* and -R* are the stellar mass and radius, 
respectively. The solar escape velocity is u^q — 617 km s~^ . 
Based on these assumptions, Table[l]shows electron number 
densities and solar mass loss rates, Mq — iiempUeQ ■ 47rr^, 
as a function of height above the solar surface near helmet 
streamers and in coronal holes, taken from |Guhathakurta| 
et al. (19961. It can be seen that all values for the solar 



mass loss rates are within the range of the observations. 
The values are higher than the value of 2 • 10~^* Mq yr~^ 
since the actual solar Alfven surface is probably higher than 
5Rq , where the density is lower. 

To derive a more general expression for stellar mass loss 
rates, we can simply equalize the squares of the Alfven and 
escape speeds at the Alfven surface to obtain: 



with M — 'i-Rr\p{rA)ue being the mass loss rate through 
the Alfven surface, assuming it is a sphere. Defining d as 
the average distance of the Alfven surface in stellar radii 
(ta = dRi,), and assuming that the magnitude of B falls like 
d~^ so that B = Bod~^, we can finally obtain an expression 
for the stellar mass loss rate as a function of the magnetic 
and gravitational stellar parameters: 



M 



BJRj 
d-^u^ 



(3) 



The definition and interpretation of Bo can be one of the 
following. 1) The weak, continuos magnetic field on the stel- 



lar surface as observed on the Sun by Schrijver et al. ( 2011 1; 
2) the weak dipole component of the stellar field; or 3) the 
"floor" value of the open magnetic flux ( Owens et al.|2008 |. 
For the solar case, any choice will set Bq to be of the order 
of 2 — 5 G. With the choice of d ranging between 2 — IO-R0, 
and solar escape velocity, we obtain a solar mass loss rate 
ranging between 10"^''' — 10"^"^ Mq yr~^ - a range that is 
in agreement with observations. 

For Sun-like stars, Bq and d are unknown. The value of 
Mi^lR-t, is about the same for most stars and does not ex- 
ceed 2.5 for B and O type stars ( Binney fc Merrifield|1998 |. 
Therefore, stellar escape velocities of Sun-like stars are lim- 



upper panels in Figure |6] show the stellar mass loss rate as 
a function of Bq and d (based on Eq. Isl, using solar es- 
cape velocity. The upper-left panel shows the mass loss rate 
for values of Bo, which are up to about order of magnitude 
larger than the solar value. It can be seen that most mass 
loss rates are in the range 10"^^ — 10~^^ Mq yr~^ . The 
rest of the plot covers unrealistic values, since d gets larger 
as Bq increases. The upper-right panel shows similar plot 
for much stronger Bo and larger value of d, which represent 
young active stars. In this plot, the mass loss rates range 
between 10"^* - 10"® Mq yr'^. 

To consider a realistic range, one should keep in mind 
that the extreme strong fields of up to few Kilo-Gausses 
observed on active stars are rather local, and are not the 
average weak component nor the dipole field. In order for the 
stellar wind to open the magnetic field lines and allow the 
mass to escape, the dynamic pressure of the wind, pd = pu^, 
must overcome the field's magnetic pressure, ps = Bo/S-Kd^. 
In the solar case, the dynamic pressure of the solar wind near 
the Alfven surface is of the order of 10~^ — 10"'* dyne cm~^ . 
The lower plot of Figure [6] shows the magnetic pressure as a 
function of Bo and d for the strong field range. It can be seen 
that the typical solar dynamic pressure of 10~* dyne cm~^ 
is not seen below d = lOi?*. Therefore, even for stronger 
fields, realistic mass loss rates are within the range of lO"'^^ — 



lO"'-" Mq yr- 
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4 CONCLUSIONS 

I study the dependency of solar mass loss rate on the level of 
solar activity and solar X-ray flux in order to test whether 
X-ray flux can be a good proxy for stellar mass loss rates. I 
exclude short-term transient and spikes from the solar wind 
and solar X-ray data, in order to better estimate the ambient 
values of the solar mass loss rate and solar X-ray flux at 
a particular period of time. I use solar wind data taken at 
different position in the solar system during different periods 
of solar activity. 

Observations show that the solar mass loss rate is scat- 
tered around the value of 2 • 10^^'* Mq yr~^ with variations 
of a factor of 2-5, while the solar X-ray flux is modified by an 
order of magnitude or even higher through the solar cycle. 
The scatter of the solar mass loss rate is wider for the slow 
wind and is more narrow for the fast wind, since the slow 
wind is more sporadic than the rather steady fast wind. The 
relative constancy of the solar mass loss rate is due to the 
fact that it is determined by the rather constant solar open 
magnetic flux. Solar activity is governed by the closed flux, 
which changes much more dramatically through the solar 
cycle. 

I derive a simple expression for stellar mass loss rates 
as a function of the stellar ambient weak magnetic fleld, the 
stellar radius, the stellar escape velocity, and the average 
height of the Alfven surface. This expression predicts stellar 



mass loss rates of 10 



10" 



Mq yr ^ for Sun-like stars. 
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Figure 1. Solar mass loss rate as a function of X-ray flux for all Ulysses data (top) and observation near Earth (bottom) taken between 
1996-2006. 
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Figure 2. Solar wind speed (left) and number density (right) as a function of X-ray for Ulysses data (top) and observation near Earth 
(bottom) taken during solar minimum (1996, blue) and solar maximum (2000, red). 
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Figure 3. Top two panels - solar mass loss rate (left) and solar kinetic energy (right) as a function of X-ray flux for Ulysses data (top) 
and observation near Earth (bottom) taken during solar minimum (1996, blue) and solar maximum (2000, red). Bottom - solar mass loss 
rate (left) and solar kinetic energy (right) as a function of count through the year taken during solar minimum (1996, blue) and solar 
maximum (2000, red). 



1965-2010 



— 1.E-13 r 



g 1.E-15 

^ 0.00 



— 1.E-13 s 
o 
g l.E-14 E^ 

O ; 

g 1.E-15_ 



5.00 10.00 15.00 

X-roy Flux [ergs s"' cm"^] 



6.E-07 


' ' ' ' 








1 




- 




- Y=-5.8- 10"'+9.8-10 


-'.> 


-5.3-10"".Xi 






X 




* ilC 

XX X 

X ^xx 
* X* X 


4.E-07 
2.E-07 




X 




-h3 k 

X * x\ 
- x'^ 


+ 


+ 


- 


O.E + 00 


dSw* , , 1 






, , * , X , 


1 







50 



200 



100 150 

Sunspot Number 

Voyoger Doily Averoged Doto 1980-2009 



50 100 150 

Doily Averoge Sunspot Number 

Voyoger Monthly Averoged Doto 1980-2009 




1995 
Yeor 



Figure 4. Top-left - long-term solar mass loss rate as a function of sunspot number and a scaled X-ray flux for data taken between 
1965-2010 (monthly averages). Top-right - the scaling function of the X-ray flux with the sunspot number used in the top-left plot. 
Bottom - solar mass loss rate as a function of time for Voyager II data averaged over a day (left) and a month (right). 
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Figure 5. Top and middle panels: solar mass loss rate as a function of wind kinetic energy (botii in logarithmic scale) for Ulysses 
(top) and 1 AU (middle) date taken during 1996 (left) and 2000 (right). Bottom-left: similar plot for date taken near 1 AU between 
1965-2010. Bottom-right: similar plot for data taken by Voyager II between 1980-2009. Blue dots are of the solar wind, while red dots 
are of measurements of the Heliosheath. 
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Figure 6. Top: logarithm of the stellar mass loss rate based on Eq. plfor smaller values of Bq and d (left) and for larger values of Bq 
and d (right). Bottom: logarithm of pB in dyne cm~^ as a function of Bq and d. 



